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We discuss the colors of 99,088 galaxies selected from the Sloan Digital Sky Survey 
Data Release 1 "main" spectroscopic sample (a flux limited sample, rp et < 17.77, 
for 1360 deg 2 ) in the rest-frame Stromgren system (uz, vz, bz, yz). This narrow-band 
(^200 A) photometric system, first designed for the determination of effective tem- 
perature, metallicity and gravity of stars, measures the continuum spectral slope of 
galaxies in the rest-frame 3200-5800 A wavelength range. We synthesize rest-frame 
Stromgren magnitudes from SDSS spectra, and find that galaxies form a remarkably 
narrow locus (~0.03 mag) in the resulting color-color diagram. Using the Bruzual & 
Chariot population synthesis models, we demonstrate that the position of a galaxy 
, along this locus is controlled by metallicity and age of the dominant stellar popula- 

ce ■ tion. The distribution of galaxies along the locus is bimodal, with the local minimum 

corresponding to an ~1 Gyr old single stellar population. The position of a galaxy per- 
pendicular to the locus is independent of metallicity and age, and reflects the galaxy's 
dust content, as implied by both the models and the statistics of IRAS detections. 

Comparison of the galaxy locus in the rest-frame Stromgren color-color diagram 
with the galaxy locus in the Hg-D n (A000) diagram, utilized by Kauffmann et al. (2003) 
to estimate stellar masses, reveals a tight correlation, although the two analyzed 
spectral ranges barely overlap. Furthermore, the rest-frame r — i color (5500-8500 A 
wavelength range) can be "predicted" with an rms of 0.05 mag using the rest-frame 
Stromgren colors. Overall, the galaxy spectral energy distribution in the entire UV 
to near-IR range can be described as a single-parameter family with an accuracy of 
0.1 mag, or better. This nearly one-dimensional distribution of galaxies in the multi- 
dimensional space of measured parameters strongly supports the conclusion of Yip et 
al. (2004), based on a principal component analysis, that SDSS galaxy spectra can be 
described by a small number of eigenspectra. Furthermore, the rest-frame Stromgren 
colors correlate tightly with the classification scheme proposed by Yip et al. (2004) 
based on the first three eigenspectra. Apparently, the contributions of stellar popula- 
tions that dominate the optical emission from galaxies are combined in a simple and 
well-defined way. We also find a remarkably tight correlation between the rest-frame 
Stromgren colors of emission-line galaxies and their position in the Baldwin-Phillips- 
Terlevich diagram. These correlations between colors and various spectroscopic diag- 
nostic parameters support earlier suggestions that rest-frame Stromgren photometry 
© 0000 RAS offers an efficient tool to study faint cluster galaxies and low surface brightness objects 

without performing time-consuming spectral observations. 
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1 INTRODUCTION 

The studies of galaxies have been recently invigorated due 
to the advent of modern sensitive large-area surveys across 
a wide wavelength range. The Sloan Digital Sky Survey 
(SDSS, York et al. 2000, Stoughton et al. 2002, Abazajian 
et al. 2003) stands out among these surveys because it has 
already provided UV to near-IR five-color imaging data, and 
high-quality spectra for over 100,000 galaxies. The spec- 
troscopic galaxy sample is defined by a simple flux limit 
(Strauss et al. 2002), and will eventually include close to 
1,000,000 galaxies. 

A number of detailed galaxy studies based on SDSS 
data have already been published. Strateva et al. (2001) 
and Shimasaku et al. (2001) demonstrated a tight corre- 
lation between the u — r color, concentration of the galaxy's 
light profile as measured by the SDSS photometric pipeline 
photo (Lupton et al. 2002), and morphology. Blanton et al. 
(2001) presented the SDSS galaxy luminosity function, and 
Kauffmann et al. (2003ab) determined and analyzed stellar 
masses and star-formation histories for 100,000 SDSS galax- 
ies. A large number of other, published and ongoing, studies 
are based on the rich dataset provided by SDSS. 

SDSS galaxy spectra have high quality (R~1800, 3800- 
9200 A wavelength range, spectrophotometric calibration 
better than 10%; for more details see Stoughton et al. 
2002), and have been used in a number of detailed inves- 
tigations, including studies of elliptical (Eisenstein et al. 
2003, Bernardi et al. 2003abcd), star- forming (Hopkins et 
al. 2003), and active galaxies (Kauffmann et al. 2003c, Hao 
et al. 2005ab, Zakamska et al. 2003, 2004). Yip et al. (2004) 
analyzed the SDSS spectra of 170,000 galaxies using princi- 
pal component analysis (Karhunen-Loeve transform). They 
demonstrated that SDSS galaxy spectra can be described 
by a small number of eigenspectra: more than 99% of the 
galaxies are found on a two-dimensional locus in the space 
spanned by the ratios of the first three eigencoefficients. 
Even the most extreme emission-line spectra can be de- 
scribed to within the measurement noise with only 8 eigen- 
spectra. An efficient single-parameter classification scheme 
is proposed, and it suggests three major galaxy classes. An- 
other remarkable result of their analysis is a strong correla- 
tion between various spectral lines, encoded in, and readily 
discernible from the eigenspectra. 

It is not yet known what the results of Yip et al. imply 
for the detailed distribution of galaxies in color space, but it 
is obvious that a high degree of structure should exist. Also, 
a strong correlation is expected between colors and various 
spectroscopic parameters, such as the strength of the Hg 
line and the 4000 A break, that were recently utilized by 
Kauffmann et al. (2003a) to estimate stellar masses and the 
the dust content for SDSS galaxies. 

The main goal of this paper is to investigate whether 
the intrinsic simplicity of galaxy spectra implied by the re- 
sults of Yip et al. can be reproduced in rest-frame color-color 
diagrams, and whether the position of a galaxy in these dia- 
grams is related to parameters discussed by Kauffmann et al. 
(2003a). We present an analysis of correlations between col- 
ors and various spectroscopic diagnostic parameters, includ- 
ing model-dependent estimates of stellar masses and dust 
content. We utilize the "rest-frame" Stromgren photomet- 
ric system that measures the continuum spectral slope of 
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Figure 1. The top panel shows renormalized (see eq. [3J filter 
transmission curves for the SDSS photometric system (ugriz), 
and for the Stromgren photometric system. The bottom panel 
emphasizes the 3800-7000 A region. The two spectra are typical 
for blue and red galaxies, and the four labeled emission lines are 
used to separate star-forming from AGN galaxies. The horizontal 
bar in the top left corner marks the wavelength region used in the 
analysis by Kauffmann et al. (2003). 



galaxies in the rest-frame 3200-5800 A wavelength range. 
In principle, any set of bandpasses could be used to syn- 
thesize the colors. An additional motivation for synthesiz- 
ing galaxy colors in the rest-frame Stromgren photometric 
system is a suggestion that it is an efficient tool for study- 
ing faint cluster galaxies and low surface-brightness objects 
without performing time-consuming spectral observations 1 
(Fiala, Rakos & Stockton 1986). The unprecedented num- 
ber of high-quality SDSS spectra can be used to analyze in 
detail the advantages and limitations of that method, and 
relate galaxy colors in the Stromgren system to parameters 
commonly measured using SDSS data. 

We describe the synthesis of rest-frame Stromgren mag- 
nitudes from SDSS spectra in Section 2. In Section 3 we an- 
alyze the distribution of galaxies in the resulting color-color 
diagram, and interpret it using the Bruzual & Chariot (1993, 
2003) population synthesis models. In Section 4 we demon- 
strate and study in detail correlations between rest-frame 
Stromgren colors and various diagnostic spectroscopic pa- 
rameters. In Section 5 we correlate the rest-frame Stromgren 
colors with the spectral eigencoefficients. In Section 6 we 
study the properties of the galaxies at sufficiently large red- 
shifts (> 0.18) to have reliable synthesized uz magnitudes. 
We discuss and summarize our results in Section 7. 



1 Observations in the rest-frame Stromgren photometric system 
use filters adjusted (in hardware) to the known cluster redshift 
(say, from the brightest cluster galaxy). 
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2 THE SYNTHESIS OF REST-FRAME 
STROMGREN PHOTOMETRY 

2.1 The Stromgren Photometric System 

The Stromgren uvby narrow-band (~200 A) filter system 
was designed for measuring the temperature, chemical com- 
position and surface gravity of stars, without resorting to 
spectroscopy (Stromgren 1966, for a recent compilation of 
over 100,000 measurements for ~63,000 stars see Hauck & 
Mermilliod 1998) . The bandpasses bracket the 4000 A break 
and cover three regions in the 3200-5800 A spectral region, 
which makes it a powerful tool for investigating stellar pop- 
ulations in star clusters, or composite systems such as galax- 
ies. 

The filter system used here (uz, vz, bz, yz; see next Sec- 
tion for details) is a somewhat modified Stromgren system, 
such that the filters are slightly narrower, and the filter re- 
sponse curves are more symmetric, than in the original sys- 
tem (Odell et al. 2002). The system characteristics are (see 
Fig.GJ: 

• The uz filter (A e // = 3500 A) is shortward from the 
Balmer discontinuity (4000 A break), and provides a mea- 
sure of hot stars due to recent star formation. 

• The vz filter (A e // = 4100 A) extends bluewards from 
4600 A, but redwards from the Balmer discontinuity. This 
wavelength region is strongly influenced by metal absorp- 
tion lines, particularly for spectral classes F and M which 
typically dominate the galaxy light. The uz — vz color is a 
good measure of the 4000 A break. 

• The bz (A e // = 4675 A) and yz (A e// = 5500 A) fil- 
ters extend redwards from 4600 A, where the influence of 
absorption lines is small. The bz — yz color is thus a good 
measure of the temperature color index, practically free of 
metallicity and surface gravity effects. For old stellar pop- 
ulations, the vz — yz and bz — yz color indices presumably 
serve as age and metallicity indicators. 



2.2 The Rest-frame Stromgren Photometric 
System 

The rest-frame Stromgren photometric system utilizes filters 
whose bandpasses are redshifted and stretched to correspond 
to the Stromgren bands in the rest-frame. This method, 
introduced by Fiala, Rakos & Stockton (1986), alleviates 
the need for K— corrections, and thus allows robust color 
comparison for galaxies at different redshifts. For example, 
Rakos & Schombert (1995) studied the color evolution of 17 
clusters of galaxies spanning the redshift range from 0.2 to 
0.9. They demonstrated that the colors of the red popula- 
tion are consistent with passive evolution models for a single 
stellar population which formed around z = 5. On the other 
hand, the fraction of blue cluster galaxies dramatically in- 
creases from 20% at z = 0.4 to 80% at z = 0.9. The rest- 
frame Stromgren photometric system also provides a good 
method for detecting cluster members, and enables an ef- 
ficient spectrophotometric classification of galaxies (Rakos, 
Dominis and Steindling 2001; Rakos & Schombert 2004; and 
references therein). 
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Figure 2. The comparison of imaging and synthesized SDSS 
magnitudes (colors) for stars, with the latter tied to imaging 
point-spread-function magnitudes. The first column displays the 
color or magnitude difference as a function of magnitude; the 
large symbols connected by lines show the ±3ct envelope. The 
distribution of color or magnitude differences at the bright end 
is displayed in the middle column. The last column shows the 
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Figure 3. Analogous to Fig. [5] except that the results are shown 
for galaxies, and the spectrophotometric magnitudes are tied to 
imaging model magnitudes. 
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2.3 The Synthetic Stromgren Photometry from 
SDSS spectra 

The synthetic Stromgren magnitudes are synthesized on the 
AB system 2 (Oke & Gunn 1983). An AB magnitude, m, in 
a bandpass S(u), is defined by 



fd(log 10 v)f v (v)S(») 

-2.51og 10 jr— — — 48.6, 

J d(log 10 v)S{v) 



(1) 



where /„ is the specific flux per unit frequency (for more 
details see Fukugita et al. 1996). The zero-point (48.6) is 
chosen such that an object with a specific flux of 3631 Jy 
has m = (i.e. an object with /„=const. has an AB magni- 
tude equal to the Johnson V magnitude at all wavelengths). 
The function S(u) includes the response of the detector, 
the transmissivity of the filter, the telescope, and the at- 
mosphere of the Earth at some representative air-mass. It 
is convenient to express the above equation in terms of the 
specific flux per unit wavelength (the form used to report 
SDSS spectra, Stoughton et al. 2002) as 



m = -2.51og 10 / <p*(\)f x (\)dA 



(2) 



where we introduced a quantity, hereafter called the 
normalized filter. 



10 1 



l XS(X) 



:J X- 1 S{X)dX 



(3) 



Computing m is, therefore, reduced to a simple integral 
of the SDSS spectrum (with the spectrum given in 10~ 17 
erg/cm 2 /s/A, and A in A) with a pre- computed normalized 
filter. The normalized filter curves are shown in Fig. 

SDSS spectra are obtained using a plate with 640 spec- 
troscopic fibers (diameter 3", corresponding to ~6 kpc at 
the redshift of 0.1), which are fed to two independently cal- 
ibrated spectrographs (Schlegel et al., in prep). Spectra are 
corrected for the Galactic interstellar dust reddening dur- 
ing calibration. As shown by Vanden Berk et al. (2002), the 
accuracy of synthesized magnitudes for point sources can 
be improved by requiring that the mean offset between the 
imaging (g, r, i) and synthesized magnitudes is zero for each 
spectrograph (this information is not available when cali- 
brating spectra). We confirmed their result for point sources, 
and found that it also improves the synthesized photome- 
try for galaxies (as inferred from the width of the galaxy 
color locus, discussed below). The synthesized photometry 
for galaxies can be further improved by using galaxies (in- 
stead of point sources) to correct for the mean offsets be- 
tween the imaging (we use model magnitudes, see Stoughton 
et al. 2002 for details) and synthesized magnitudes. Here- 
after, we adopt the latter method. When correcting synthe- 
sized Stromgren magnitudes, we linearly interpolate (as a 
function of wavelength) the g and r band corrections. 



2.4 The Tests of Synthesized Magnitudes 

We synthesized SDSS magnitudes {g, r, i) and rest-frame 
Stromgren magnitudes (uz, vz, bz, yz) 3 for 99,088 unique 
"main" (rp et < 17.77) galaxies from SDSS Data Release 1 
(Abazajian et al. 2003), as well as for stars from the same 
spectroscopic plates. The accuracy of synthesized magni- 
tudes (relative to imaging magnitudes, which are accurate 
to 0.02 mag, Ivezic et al. 2003) for a subsample of stars is 
summarized in Fig. |2] The colors are reproduced with the 
expected accuracy (~0.08 mag), while there are systematic 
errors in the overall magnitude scale at the level of 0.10 mag 
(root-mean-square scatter). These results are in agreement 
with Vanden Berk et al. (2004). 

Fig. |3] summarizes the obtained accuracy for galaxies. 
Again, the colors are reproduced with the expected accu- 
racy, while systematic errors in the overall magnitude scale 
are ~0.2 mag. Note that using galaxies (instead of point 
sources) to correct for the offset between imaging and syn- 
thesized magnitudes introduces a systematic trend in the 
Ar distribution at the faint magnitude end. Nonetheless, as 
evident from the top and bottom panels, this does not signif- 
icantly affect the synthetic colors as the overall trend nearly 
cancels out. In addition, the distribution of color-differences 
normalized by the expected errors (last column in Fig. [3] 
top and bottom panels) shows that the errors are estimated 
correctly (the width of the distributions is ~1). Therefore, 
we do not expect this to affect the synthetic photometry 
significantly. The errors in rest-frame Stromgren colors are 
expected to be smaller than in g — r and r — i colors be- 
cause the relevant wavelength range is shorter (and thus less 
sensitive to large-scale errors in the spectrophotometric cal- 
ibration of SDSS spectra). We estimated errors in vz — yz 
and bz — yz colors by comparing measurements for 1253 
galaxies observed twice: the scatter (root-mean-square de- 
termined from the interquartile range) for the vz — yz color 
is 0.09 mag, and for bz — yz color 0.05 mag. We tested for 
plate-to-plate systematics by analyzing the position of the 
red peak in the color distribution of galaxies (see the next 
Section). This position shows a scatter of ~ 0.02 mag for 
vz — yz color, and ~ 0.01 mag for bz — yz color. We con- 
clude that the synthesized colors are sufficiently accurate 
for a detailed investigation of rest-frame colors. This is the 
first time that rest-frame optical colors with an accuracy of 
<0.1 mag are available for a sample of ~ 100,000 galaxies. 



3 THE COLORS OF GALAXIES IN THE 
REST-FRAME STROMGREN SYSTEM 

In this Section we analyze the distribution of galaxies in 
the rest-frame Stromgren color-color diagram, discuss pos- 
sible systematic errors and interpret the colors using model 
spectra. 



2 Note that all the published work by Rakos and collaborators, 
who pioneered the redshifted Stromgren photometric system, is 
expressed on the Vega system. 



3 We discuss the implications of the uz magnitude separately in 
Section [H] since it can only be synthesized for galaxies with red- 
shifts beyond ~0.18. 
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for the red side, with ci = vz — yz and c% — bz — yz. 
The coefficients are normalized such that the errors in PI 
and P2 are the same as the errors in vz — yz and bz — yz, in 
the limit when the latter two are the same. It is noteworthy 
that the P2 axis is nearly parallel to the color index mi = 
(vz — yz) — 2 (bz — yz), used as a metallicity estimator for 
single stars (e.g. Twarog 1980). 

The distributions of PI and P2 colors are shown in 
Fig. |S] by thick lines. The width of the locus in the P2 di- 
rection is surprisingly small 4 (0.03 mag, and 0.32 mag in 
the PI direction). Indeed, the rest-frame colors for a large 
number of galaxies have not been measured to date with 
such an accuracy. This small value both testifies that the 
errors in synthesized colors are small, and that the slope 
of galaxy spectra in the 4000-5800 A wavelength range is a 
nearly one-parameter family. 



Figure 4. The distribution of 99,088 SDSS "main" galaxies in the 
rest-frame Stromgren color-color diagram. Each galaxy is repre- 
sented by a dot, and the overall distribution is outlined by linearly 
spaced isopleths. The principal axes PI and P2 are defined for the 
two regions separated by vz — yz = 0.646. 
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Figure 5. The color probability distributions of the principal 
rest-frame Stromgren colors for SDSS "main" galaxies (note dif- 
ferent scales for x axis). The thick lines show the distributions 
for the whole sample, and thin lines for the subsample of galaxies 
detected by IRAS. The a values shown in each panel are root- 
mean-square scatter (determined from the interquartile range) 
for the whole sample. 

3.1 The Galaxy Locus and Principal Axes 

The distribution of 99,088 galaxies in the rest-frame 
Stromgren color-color diagram is shown in Fig.|I] The distri- 
bution resembles a remarkably narrow locus. To quantify its 
width, we define a set of principal axes (PI, P2) where PI 
measures the position along the locus, and P2 perpendicular 
to it. Since the locus is not perfectly straight, we separately 
fit the blue (vz — yz < 0.6) and red (vz — yz > 0.7) ends. 
As the boundary between the two (P1,P2) definitions we 
choose vz — yz = 0.646, where the two fits intersect. We 
obtain for the blue side 



PI = 0.911 (ci - 0.646) + 0.412 (c 2 - 0.261) 
P2 = -0.412 (ci - 0.646) + 0.911 (c 2 - 0.261), 
and 

PI = 0.952 (ci - 0.646) + 0.307 (c 2 - 0.261) 
P2 = -0.307 (ci - 0.646) + 0.952 (c 2 - 0.261), 



(4) 
(5) 

(6) 
(7) 



3.2 Fiber Aperture Effects on PI & P2 

SDSS spectra are obtained with a fixed 3 arcsec diameter 
fiber positioned as close as possible to the center of the 
galaxy. Since the fiber aperture may not include the entire 
galaxy, it is possible that colors derived from SDSS spectra 
are biased red by the bulge contribution. This bias should 
be the largest for intermediate- and late-type galaxies where 
the fiber aperture may miss a significant fraction of the outer 
disk and its star formation. While this effect was discussed 
in detail elsewhere (e.g. Kauffmann et al. 2003a and Kewley 
et al. 2005), here we explore its impact on the PI and P2 
colors. 

The fiber aperture effect will tend to make the PI color 
redder for galaxies with large angular size because bulges 
are generally redder than disks. The expected effect for the 
P2 color is not so obvious. We argue below (see Section |3~H 
that the P2 color is related to dust, and if this is true, the ex- 
pected effect is for the P2 color to become bluer for galaxies 
with large angular size. 

Figures [U and [7| show the dependence of the PI and P2 
colors on the Petrosian 50% radius in the z band 5 for galax- 
ies with the r band absolute magnitudes in the range -19 to 
-18. This magnitude slice is selected because it is dominated 
by blue galaxies where the expected effect is the strongest 
(see Fig. I2H . As the Petrosian radius increases, the bulge 
contribution to the flux captured by the fiber also increases. 
As evident, small but measurable dependence of the PI color 
on the Petrosian radius is indeed present. Given the observed 
range of the Petrosian radius, the expected bias in the PI 
color is <0.05 mag, with the root-mean-square scatter con- 
tribution to the PI measurement error of ~0.02 mag. Com- 
pared to the observed PI range, this is a small effect. The 
fiber aperture effect on the P2 color is entirely negligible. 
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Figure 6. A test for aperture effects, due to 3 arcsec SDSS spec- 
troscopic fiber diameter, on the PI color. The linearly spaced 
contours show the distribution of galaxies with absolute Pet- 
rosian magnitude in the r band between -19 and -18 in the PI 
color vs. Petrosian 50% radius in the z band (an estimator of 
the galaxy's size) diagram. The aperture effect is expected to 
be the strongest for blue galaxies with PKO. The large circles 
show median PI colors for bins of the Petrosian radius in the z 
band, computed for galaxies with PKO. The middle line is a best 
straight line fit to these medians, and the other two lines illus- 
trate 2<r deviation around the median. The slope of the best fit 
line is 0.015 mag/arcsec. 
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Figure 8. The Bruzual & Chariot model predictions for the de- 
pendence of the principal rest-frame Stromgren colors on the age 
and metallicity ([Z]=0.004, squares; 0.02, circles; 0.05, dots) of 
single stellar populations. Note the different scales for the two 
colors. The histograms on the y axes show the observed color 
distributions. 
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Figure 7. Analogous to Fig. [H] except for the P2 color and 
only for galaxies with PKO. The slope of the best fit line is - 
0.0005 mag/arcsec. 



3.3 The Stellar Population Model Colors 

The principal axes (P1,P2) are defined by the observed mor- 
phology of the distribution of galaxies in the rest-frame 
Stromgren color-color diagram. Here we relate these axes to 
the stellar age and metallicity using the Bruzual & Chariot 
(1993, 2003) population synthesis models. 

Using the same method as for observed spectra, we 
synthesize rest-frame Stromgren colors from model spectra. 
Fig. |H| shows the model colors as a function of the age of a 
single stellar population, and for three representative met- 
alicities. A single stellar population corresponds to a single 
epoch of star formation, while real galaxies likely have more 
complex star formation histories. In addition, for single stel- 
lar population models the star-formation rate history, which 
is an important discriminating parameter (e.g. Kennicutt 
1998), is unavailable. Nevertheless, the color range spanned 
by a single stellar population provides a robust constraint on 
the possible color distribution for any star formation history 
(because it brackets the possible color range). 

The observed range of the PI color can be explained by 
stellar populations with ages from 100 Myr to 10 Gyr. Note 
that PI > corresponds to populations older than ~1 Gyr, 



4 It may be surprising that the width of the P2 color distribution 
is smaller than the errors of the colors used to synthesize it. This 
is due to covariances in the vz — yz and bz — yz errors. 

5 For more details about the Petrosian radii see Stoughton ct al. 
2002, and Strauss et al. 2002. 
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but the age estimate is very uncertain due to the unknown 
metallicity. 

Models provide a simple explanation for the observed 
small locus width: the P2 color is a color projection which is 
nearly independent of both age and metallicity 6 . At any age, 
P2~0 to within 0.02-0.03 mag, and for populations older 
than ~1 Gyr, P2~0.01 to within 0.01 mag. That is, what- 
ever is age or metallicity, P2 is always to within a few 
hundreths of a magnitude. This is very different from PI 
which varies by 2 magnitudes as a function of age, and up 
to 0.6 mag as a function of metallicity. The ability of models 
to explain the mean value and the very narrow distribution 
of the P2 color is a considerable success because the princi- 
pal axes PI and P2 are defined solely by the morphology of 
the data distribution in the bz — yz vs. vz — yz color-color 
diagram. 

The observed width of the P2 distribution is wider 
(0.03 mag) than the model prediction (~0.01 mag). While 
this difference could be easily attributed to the measure- 
ment scatter or to the simplicity of the models, we find that 
the P2 color is correlated with the estimate of the galaxy 
dust content, A z (the effective dust extinction in the SDSS z 
band), determined by Kauffmann et al. (2003a) and further 
discussed in Section [4.21 Hence, the larger scatter observed 
in the P2 color than predicted by the model is, at least par- 
tially, due to dust reddening (which is not included in the 
models) . 

It could be argued that the P2 — A z correlation is simply 
due to the fact that both quantities are determined from 
SDSS spectra. However, they are determined using different, 
barely overlapping, spectral ranges. While Kauffmann et al. 
analysis utilized the 3850-4160 A spectral range to constrain 
their models (and then compared imaging g — r and r — i 
colors to model-predicted colors to estimate dust reddening) , 
P2 is determined from the 4000-5800 A wavelength range 
(see the bottom panel in Fig.0. 

Another reasonable objection is that the A z estimates 
are, of course, model dependent and consequently may not 
be a good measure of the dust content. However, Obric et 
al. (2006) demonstrate that galaxies detected by IRAS have 
systematically higher values of A z . This finding strongly sug- 
gests that A z does provide an estimate of the dust content. 
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Figure 9. The comparison of P2 distributions for all galaxies and 
for the subsample detected by IRAS, for 0.2 mag wide bins in PI, 
in the range -0.7 to 0.7 (row by row, from top left to bottom right). 
The thin lines show the distributions for the whole sample, and 
the thick lines for the IRAS-detected subsample. The A values 
shown in each panel are the difference between the medians for 
the two distributions, and are also shown as a function of PI in 
the bottom right panel. 



3.4 The Colors of IRAS-detected galaxies 

Since the P2 — A z correlation implies that P2 is also a mea- 
sure of the galaxy dust content, galaxies detected by IRAS 
should have systematically higher values of P2. Using the 
same SDSS-IRAS galaxy sample discussed by Obric et al. 
(~2200 probable matches for the sample of 99,088 galaxies 
discussed here), we compared the PI and P2 distributions 
for IRAS-detected galaxies to those for the full sample. As 
evident from Fig. [B] galaxies detected by IRAS have system- 
atically higher P2 values. 

6 The P2 color, while much less sensitive to metallicity than the 
PI color, does become redder with increasing metallicity for pop- 
ulations older than 1 Gyr. However, the effect is exceedingly small 
(0.03 mag per dex). For individual stars (P2 is nearly parallel to 
mi, c.f. section lij.ll this is large enough because their effective 
temperature is typically constrained much better than the star- 
formation history of a galaxy. 



The linear definition of the principal axes PI and P2 
(eqs. 0D~II1) is not a perfect description of the galaxy locus 
shown in Fig. 2] because the locus is curved. Hence, the P2 
distribution slightly changes along the locus. Since IRAS- 
detected galaxies are biased towards blue galaxies, it could 
be that the difference in P2 distributions for IRAS-detected 
and all SDSS galaxies shown in Fig.|5]is caused by the vari- 
ation of the P2 distribution along the locus. Fig. [5] demon- 
strates that this is not the case: for each segment of the 
locus defined by a narrow bin of PI, IRAS-detected galax- 
ies have systematically higher P2 colors. This offset strongly 
suggests that P2 is a measure of the galaxy dust content. 

The shift of the P2 distributions between SDSS-IRAS 
and all SDSS galaxies is about the same as the P2 distri- 
bution width, so P2 is a rather noisy measurement of the 
dust content. Nevertheless, we demonstrate below that it is 
sufficiently efficient to study differences in the dust content 
between various galaxy populations separated using emis- 
sion lines. 
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4 THE CORRELATIONS BETWEEN THE 

REST-FRAME STROMGREN COLORS AND 
SPECTRAL PARAMETERS 

In the previous Section we have shown the connection be- 
tween the two principal rest-frame Stromgren colors and 
galaxy physical parameters: PI depends, in a degenerate 
way, on stellar age and metallicity, while P2 is correlated 
with the dust content, and independent of age and metallic- 
ity. In this Section we further analyze correlations between 
Stromgren colors and various spectral parameters. We show 
that the position of an emission line galaxy in a diagram 
commonly used to separate star-forming from AGN galaxies 
is well correlated with the PI color, that PI and P2 colors 
are correlated with the parameters discussed by Kauffmann 
et al. (2003ab), and that the rest-frame r — i color can be pre- 
dicted to better than 0.05 mag from the PI color. All these 
correlations suggest that galaxy spectra are simple and well- 
defined superpositions of stellar spectra. 



4.1 The Baldwin-Phillips- Terlevich Diagram 

The BPT diagram (Baldwin, Phillips and Terlevich, 1981) 
is a standard method for classifying emission-line galaxies as 
star-forming and AGN galaxies. With the advent of SDSS 
data, it is now possible to construct the BPT diagram for 
an unprecedented number of galaxies (e.g. Ivezic et al. 2002, 
Kauffmann et al. 2003c, Hao et al. 2005ab, Heckman et al. 
2004). Obric et al. (2006) have shown that the presence of 
emission lines in a galaxy's spectrum is well correlated with 
the u — r color and the concentration index determined solely 
from the imaging data. Galaxies without emission lines tend 
to have larger concentration indices and redder u — r colors 
than galaxies with emission lines. Furthermore, the distri- 
bution of emission-line galaxies in the BPT diagram is also 
correlated with the u — r color and concentration index. 
Galaxies classified as star-forming have predominantly blue 
u — r colors and small concentration indices, while AGN 
galaxies have redder u — r colors and large concentration 
indices. 

Here we investigate the correlation between the strength 
of the emission lines used to construct the BPT diagram, 
and rest-frame Stromgren colors. To classify a galaxy as an 
emission-line galaxy, we require a 3<r significant detection of 
H a , Hp, [Nil] and [O77/5007] lines (the line strengths are 
determined as described by Kauffmann et al. 2003c). The 
bottom panel in Fig. 1101 shows the distribution of 43,281 
emission-line galaxies, from the sample discussed here, in 
the BPT diagram. To visualize the correlation with PI and 
P2 colors, the dots are two-dimensionally color coded ac- 
cording to their position in the P1-P2 color-color diagram 
shown in the top panel. There is a strong correspondence 
between the position of a galaxy in the BPT diagram and its 
PI color. Galaxies in the "star-forming branch" with small 
[NII]/H a ratios, for a given [Oil 1 5007]/ Hp ratio, have pre- 
dominantly blue PI colors, while AGN galaxies have redder 
PI colors. This correlation is further illustrated in Fig. 1111 
The correlation with the P2 color is not discernible. 

In the subsequent analysis, we separate emission-line 
galaxies in three groups according the their position in the 
BPT diagram: AGNs, star-forming, and "unknown". The 
adopted separation boundaries are shown by the dashed 




-1.5 -1 -0.5 0.5 

[Nll]/H a 

Figure 10. The top panel shows the distribution of SDSS "main" 
galaxies in a diagram constructed with the principal rest-frame 
Stromgren colors. Each dot corresponds to one galaxy; their dis- 
tribution is outlined by linearly spaced isopleths. The color code 
is determined by the position in this diagram, and is used in 
subsequent figures to visualize correlations of other quantities 
with PI and P2. The bottom panel shows the distribution of 
galaxies in the Baldwin-Phillips- Terlevich diagram, constructed 
with emission-line strength ratios. The dashed lines separate the 
regions populated by star-forming and AGN galaxies. They are 
motivated by the modeling results from Kewley et al. (2005) and 
are fine-tuned to isolate the peak of the observed distribution. 
For galaxies in the wedge-shaped regions the classification is un- 
certain. The dots are colored using PI and P2; note the strong 
correlation between the PI color and [Nil]/ Ha ratio, although 
these parameters are determined in different wavelength ranges, 
and measure the continuum slope and line strengths, respectively. 

lines. The last category is found at the joint of the two 
branches, and it is not obvious from the displayed data to 
which class these galaxies belong. Fig. 1121 compares the PI 
distributions of these classes. The majority of galaxies with- 
out emission lines and AGN galaxies have PI > 0, while 
star-forming galaxies have PI < 0. More than ~2/3 of galax- 
ies with PI > 0.3 are galaxies without emission lines, and 
more than ~2/3 of galaxies with PI < —0.4 are star-forming 
galaxies. At least 20% of the galaxies with PI > harbor 
an AGN. The fraction of AGN galaxies drops to essentially 
zero by PI ~ —0.3 (or by PI ~ —0.5 if all unclassified 
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Figure 11. The correlation between the first principal axis in 
the rest-frame Stromgren color-color diagram, and the [Nil]/ Ha 
line strength ratio. The distribution of all emission-line galaxies 
in the sample is shown by linearly spaced isoplcths. 
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Figure 12. The distributions of PI, the first principal axis in 
the rest-frame Stromgren color-color diagram, for subsamples of 
galaxies separated using emission line strengths: galaxies with- 
out emission lines (NOEML), star- forming galaxies (SF), AGNs, 
and emission-line galaxies which cannot be reliably classified (un- 
known). The top panel displays counts, and the bottom panel 
displays the fraction of each subsample in the whole sample. 



galaxies are AGNs) . This dramatic absence of AGNs among 
the bluest galaxies appears not to be a selection effect, as 
discussed by Kauffmann et al. (2003c, Section 4.1). 

The comparison of P2 distributions for the three sub- 
classes of galaxies is shown in Fig. 1131 As expected, AGNs 




Figure 13. The distributions of P2, the second principal axis in 
the rest-frame Stromgren color-color diagram, for subsamples of 
galaxies separated using emission line strengths: galaxies without 
emission lines (NOEML), star-forming galaxies (SF) and AGNs. 
The top panel displays counts, and the bottom panel displays 
the fraction of each subsample in the whole sample (emission-line 
galaxies which cannot be reliably classified make up the missing 
fraction) . 



and star-forming galaxies have redder P2 distributions than 
galaxies without emission lines. 

4.2 The Hg - D„(4000) Locus 

Kauffmann et al. (2003a) used the distribution of galaxies 
in the plane spanned by the strength of the Hs line and 
the 4000 A break (D„(4000)) to obtain model-dependent 
estimates of stellar masses and the dust content for SDSS 
galaxies. Given the position of a galaxy in the _r7^-D n (4000) 
plane, the most probable mass-to-light ratio is drawn from a 
model library. With the measured luminosity, this ratio then 
yields the stellar mass. The observed luminosity is corrected 
for the dust extinction determined by comparing observed 
imaging g — r and r — i colors to model-predicted colors 
(the latter do not include the effects of dust reddening). 
The reddening correction needed to make models agree with 
data is interpreted as an effective optical depth, A x , due to 
the galaxy's interstellar dust. As discussed in Section 13.31 
Obric et al. (2006) demonstrate that galaxies detected by 
IRAS have systematically higher values of A z , supporting 
the notion that A z provides an estimate of the dust content. 

The analysis by Kauffmann et al. is based on the 3850- 
4160 A spectral range, marked by the horizontal bar in the 
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Figure 14. The top panel shows the distribution of SDSS galaxies 
in the plane spanned by the strength of the Hg line and the 4000 A 
break (P„(4000)). Each galaxy is represented by a point color- 
coded according to its position in the P1-P2 color-color diagram, 
shown in the top panel in Fig. 1101 Note the strong color gradient 
indicating a correlation between the position of a galaxy along this 
locus and its PI color. The overall, strongly bimodal, distribution 
is outlined by linearly spaced isopleths. The bottom panel shows 
the distribution of the same sample in the dust content (A z ) 
vs. stellar mass (M„) diagram. A z and M, are model-dependent 
estimates derived from the measured Hg, D n (4000) and imaging 
g— r and r—i colors. The color gradients indicate good correlations 
between the positions in the A z -M* and PI-P2 planes. 



bottom panel in Fig. while the rest-frame Stromgren col- 
ors discussed here are determined using the 4000-5800 A 
range. That is, the two studies are based on practically inde- 
pendent spectral ranges (the vz filter does have overlap with 
the Hg line). Here we demonstrate that despite their dif- 
ferent wavelength range, and different techniques (color vs. 
spectral line analysis), the derived inferences about galaxies 
are similar. 

The top panel in Fig. ll4l shows the distribution of galax- 
ies in the Hg-D n (4000) plane, with the dots colored accord- 
ing to the galaxy's position in the P1-P2 color-color diagram 
(see the top panel in Fig. llUI . As discussed by Kauffmann et 
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Figure 15. Analogous to Fig. 1141 except that only galaxies with- 
out emission lines are shown (contours correspond to the full sam- 
pie). 



al. (2003a), galaxies form a well-defined locus in this plane. 
The strong color gradient along the locus indicates a corre- 
lation with the PI color. The bottom panel in Fig. ll4l shows 
the distribution of galaxies in the plane spanned by the de- 
rived parameters: the dust content and the stellar mass. The 
color gradients indicate a good correlation between the posi- 
tions in the A z -M* and P1-P2 planes. In particular, galax- 
ies with small masses (log(M*)<10) are dominated by blue 
galaxies with PI < 0, while a correlation between A z and P2 
is discernible (gradient from purple to red as A z increases) 
for high-mass galaxies (log(M*)~ll). Since emission line 
properties, discussed in the previous Section, correlate with 
PI and P2 colors, it is expected that galaxies without emis- 
sion lines, star-forming, and AGN galaxies display different 
distributions in these two diagrams. This is demonstrated 
in Figs. 1151 1161 and 1 171 For example, the detailed morphol- 
ogy of the galaxy distribution shown in the bottom panel 
in Fig. 1141 can be understood as a superposition of three 
dominant populations with distinctive positions, as seen in 
the bottom panels in Figs. 1151 1161 and 1171 It is noteworthy 
that the A z -M t correlation for star-forming galaxies shown 
in the bottom panel in Fig. 1161 is consistent with a constant 
dust-to-stellar mass ratio for all galaxies. 

In order to quantify the correlations between the po- 
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Figure 16. Analogous to Fig. 1141 except that only star-formi: 
galaxies are shown (contours correspond to the full sample). 



Figure 17. Analogous to Fig. 1141 except that only AGN galaxies 
are shown (contours correspond to the full sample). 



sitions in the Hs-D„ (4000) and P1-P2 planes, we define a 
parameter rj which measures the position along the locus in 
the (4000) plane 

rj = 0.0995 ■ (D„(4000) - 1.75) - 0.995 • H8 A (8) 

Note that »y=0 corresponds to H8a=0. Fig. H8l compares the 
•q-Pl and A z -P2 distributions for galaxies with and without 
emission lines. As evident, there is good correspondence be- 
tween the two pairs of parameters, although they are mea- 
sured using barely overlapping spectral ranges. In particular, 
note the good correlation between A z and P2 for emission 
line galaxies. Analogous diagrams comparing star-forming 
and AGN galaxies are shown in Fig. 1191 It appears that the 
A z -P2 correlation is stronger for AGN than for star-forming 
galaxies. This could be explained as a consequence of the 
more homogeneous AGN population because star-forming 
galaxies span a broad range in stellar mass, and hence in 
metallicity and surface mass density. It also seems that the 
slope of the A z -P2 correlation is different for the two types 
of galaxies, but the data are too noisy to derive a robust 
conclusion. 



4.3 The Correlations of A z and P2 with the 
H a /H(j Ratio 

An often used estimator of the galaxy dust content is the 
H a / Hp ratio. Therefore, at least some degree of correlation 
should exist between this parameter and the A z and P2 
parameters. Fig. 1201 shows the dependence of log(H a / Hp) 
on P2 and A z for star-forming and AGN galaxies. Simi- 
larly to A z vs. P2 dependence shown in Fig. 1191 it appears 
that different weak correlations exist for these two types of 
galaxies. The difference in slopes, similar to the difference 
in slopes for the A z vs. P2 correlation, is possibly due to 
different dust geometries, or perhaps different dust proper- 
ties, between AGN and star-forming galaxies (note that this 
difference also exists in the Ha/Hp vs. A z diagram, that is, 
even when P2 is not used) . A detailed analysis and interpre- 
tation of these correlations is beyond the scope of this work. 
We conclude that the similarity of the top and bottom pan- 
els in Fig. 1201 suggests that P2 is a good, though somewhat 
noisier, proxy for A z . 
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Figure 18. The top panel shows the correlation between the po- 
sition of a galaxy on the locus in the Hg-D n (A000) plane (rj) and 
the PI color. The thin contours correspond to galaxies without 
emission lines, and the thick contours to emission-line galaxies. 
The bottom panel shows the A z -P2 correlation. 



4.4 The Color — Absolute Magnitude 
Distribution 

The color distribution of galaxies depends on their abso- 
lute magnitude: the fraction of red galaxies increases with 
luminosity. Furthermore, the median color of red galaxies 
becomes even redder with luminosity (Baum 1959, Faber 
1973). These long-known results are spectacularly confirmed 
by the accurate and voluminous SDSS data (e.g. Blanton et 
al. 2003, Baldry et al. 2004, and references therein). Hence, 
the distribution of the PI principal color introduced here is 
expected to be highly correlated with luminosity. 

The top panel in Fig. l21l shows the distribution of 1,704 
galaxies selected from a narrow redshift range in the abso- 
lute magnitude vs. PI color plane. We use the Petrosian r 
band magnitudes without K-correction and WMAP cosmol- 
ogy (Spergel et al. 2003) to compute M r p e t (for a descrip- 
tion of Petrosian magnitudes see Stoughton et al. 2002, and 
Strauss et al. 2002) . The strong color-magnitude relation for 
red galaxies is evident, as well as the dramatic increase in 
the fraction of red galaxies for M r p e t < —20. Indeed, red 
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Figure 19. Analogous to Fig. 1181 except that star-forming (thick 
contours) and AGN (thin contours) galaxies are shown. 



galaxies with PI ~ 0.5 show a well-defined local maximum 
in their luminosity function (at M r p et ~ —20.7). 

This figure suggests that PI = 0.25 provides a good 
separation between blue and red galaxies (u — r = 2.22 
separator from Strateva et al. 2001 corresponds to median 
PI = 0.18). The separation of less-luminous blue galaxies 
and more-luminous red galaxies is much more pronounced 
in this diagram based on the PI color, than when using other 
colors such as K-corrected broad-band SDSS colors (see e.g. 
Fig. 7 in Blanton et al. 2003). This is presumably because 
PI is a principal axis in a rest-frame color-color diagram 
constructed with narrow-band filters (the galaxy locus in 
the r — i vs. g — r rest-frame color-color diagram, discussed 
in Section [4.71 is twice as broad as the galaxy locus in the 
bz — yz vs. vz — yz plane). 

Since the PI color is correlated with numerous other pa- 
rameters, the subsamples separated using those parameters 
should have distinct distributions in the color - absolute 
magnitude diagram. As an example, we choose the sepa- 
ration into star-forming, AGN and galaxies without emis- 
sion lines. The bottom three panels in Fig. 1211 compare the 
color-absolute magnitude distribution for each subsample 
to that for the full sample. As evident, AGN and galaxies 
without emission lines are predominantly red and luminous, 



© 0000 RAS, MNRAS 000, 000-000 



The Rest-frame Colors of SDSS Galaxies 13 



0.028 < redshift < 0.032 




-0.08 



-0.04 





P2 



0.04 



0.08 




-0.5 



0.5 
A. 



1.5 



Figure 20. The comparison of \og(H a /Hn) vs. P2 (top) and 
log(H a I Ha) vs. A z (bottom) correlations for star-forming (thick 
contours) and AGN (thin contours) galaxies 



while star-forming galaxies are blue and faint. The faint blue 
galaxies without emission lines, visible in the upper middle 
panel, are galaxies with weak emission lines that did not 
fulfill the condition of the 3cr detection for all four relevant 
emission lines. 



4.5 The Color — Mass Distribution 

Kauffmann et al. (2003a) estimated stellar masses for SDSS 
galaxies by multiplying observed luminosities (corrected for 
extinction using A z ) with the stellar mass-to-light ratio in- 
ferred from the position of a galaxy in the _ff,5-D„ (4000) 
plane. The correlation of the mass-to-light ratio and param- 
eter rj, which measures the position of a galaxy along the 
Hs-D n (4000) locus, is shown in the top panel in Fig. 1221 
The dashed line is a best fit 

log(M/L z ) = 0.000630r? 3 -0.00174r7 2 +0.01396^+0.231, (9) 

which "predicts" log (M/L z ) with an rms of 0.13 (the overall 
distribution of log (M/L z ) has an rms of 0.17, i.e. a factor of 
1.5). The upper middle panel shows the mass-to-light ratio 
as a function of the PI color. The best fit, 
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Figure 21. The color-absolute magnitude diagram for SDSS 
main galaxies. The top panel shows the absolute r band mag- 
nitude as a function of the PI rest-frame Stromgren color for 
galaxies with 0.028<redshift<0.032. The remaining panels com- 
pare the distribution of all galaxies (contours) to the distribu- 
tions of three subsamples (symbols) selected using emission lines: 
galaxies without emission lines in the upper middle panel, AGN 
galaxies in the lower middle panel, and star-forming galaxies in 
the bottom panel. 
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Figure 22. The top panel shows the correlation between the 
mass-to-light ratio (in solar units) and the position of a galaxy 
on the Hg — _D n (4000) locus, where the model-dependent mass-to- 
light ratio is taken from Kauffmann et al. (2003). The dashed line 
is a best-fit third order polynomial. The upper middle panel shows 
the correlation between the mass-to-light ratio and the position 
of a galaxy on the locus in the rest-frame Stromgren color-color 
diagram. The dashed line is a best-fit third order polynomial. 
The lower middle panel shows the mass - color diagram using the 
best-fit mass-to-light ratio determined from the PI color, and 
the bottom panel shows analogous distribution using the mass- 
to-light ratio from Kauffmann et al. (2003). Note the similarity 
of the two distributions. 



log(M/L 2 ) = 0.273 PI 3 -0.328 P1 2 +0.516P1 + 0.091, (10) 

reproduces log (M/L z ) with an rms of 0.10, equal to the me- 
dian error in estimating log (M/L z ) using the Kauffmann et 
al. (2003a) method. Hence, it is possible to use PI to convert 
the observed luminosities into stellar masses with practically 
the same accuracy as the Hs-D n (4000) method employed by 
Kauffmann et al. (2003a). The bottom two panels compare 
the mass-color distributions when the mass is determined us- 
ing PI, and when the mass is taken from Kauffmann et al. 
(2003a); their similarity confirms this conclusion. We note 
that the overall appearance of the mass-color distribution, 
in particular, the strong correlation of mass and color, es- 
sentially reflects the distribution of galaxies in the color- 
luminosity diagram (because the mass-to-light ratio has a 
very small dynamic range). 

Since PI and luminosity, L z , are correlated, as de- 
scribed in the previous Section, and M/L z is a function 
of PI, M/L z and L z should also be correlated to at least 
some extent. For example, Bernardi et al. (2003b) deter- 
mined that for elliptical galaxies M/L oc L ' 14 . We find that 
a best power-law fit for all galaxies is M/L oc L 0A . It repro- 
duces log (M/L z ) with an rms of 0.20. Since this scatter is 
larger than when determining M/L from PI, we speculate 
that galaxy colors are a better indicator of the mass-to- light 
ratio, and that its correlation with galaxy luminosity is prob- 
ably a consequence of the color-magnitude relation. 

4.6 The Color — Velocity Dispersion Distribution 

Stellar velocity dispersion is automatically measured for 
each SDSS galaxy by fitting the observed spectrum with a 
linear combination of galaxy template spectra broadened by 
a Gaussian kernel (Schlegel et al., in prep.). Detailed analy- 
sis of velocity dispersions for elliptical galaxies, as measured 
by SDSS, is discussed by Bernardi et al. (2003a) and Sheth 
et al. (2003). 

We find a good correlation between the velocity disper- 
sion and the PI color, shown in the top panel in Fig. 1231 
The correlation is especially strong for velocity dispersions 
larger than the instrumental resolution of the SDSS spectra 
(60-70 km/s, log(?;disp/100 km/s)~-0.2). A similar correla- 
tion between the velocity dispersion and the g — r color for 
elliptical galaxies is discussed by Bernardi et al. (2003d). 

The correlation between the velocity dispersion and the 
PI color appears to be a consequence of the good correlation 
between the stellar mass and the velocity dispersion, shown 
in the bottom panel in Fig. 1231 We find that the function 

log (M/M ) = 10.7 + 2 log ( I ^£ 7i ). (ID 

provides a good description for both blue and red galaxies. 
For a constant M, in the range 10.5 < log (M/Mq) < 12, the 
rms scatter in the velocity dispersion is 40 km/s, indepen- 
dent of M. For a constant velocity dispersion, the rms scatter 
in log (M/Mq) decreases from 0.33 for «di sp = 100 km/s to 
0.18 for u dis p = 250 km/s. 

The similar log (M/Mq) vs. i>dis P behavior for blue and 
red galaxies can be considered as circumstantial evidence 
that the model-dependent mass-to-light ratios must be cor- 
rect to some extent. The reason is that blue (PI < 0.25) and 
red galaxies (PI > 0.25) have different observed luminosi- 
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Figure 23. The top panel shows the correlation between the PI 
rest-frame Stromgren color and the velocity dispersion. The bot- 
tom panel shows the correlation between the galaxy's stellar mass 
and the velocity dispersion, separately for galaxies with Pl<0.25 
(solid lines) and Pl>0.25 (dashed lines). The diagonal straight 
line, log(M/A4" )=lO.7+21og(i) dis p/lOO km/s), is added for illus- 
tration 



ties (median z band luminosity for red galaxies is larger by 
a factor of 2.6 than for blue galaxies), and different velocity 
dispersions, and only the "correct" mass-to-light ratios (me- 
dians are 1.2 and 1.9 for blue and red galaxies, respectively) 
will place the inferred stellar masses on the same M — Wdisp 
relation. If the same mass-to-light ratio is used for all galax- 
ies, the inferred mass for red galaxies would be smaller by 
a factor of 1.6 than for blue galaxies with the same veloc- 
ity dispersion. That is, the two distributions shown in the 
bottom panel in Fig. 1231 would be vertically offset by ~0.2. 
The data constrain the offset between the two distributions 
to <0.1. 



4.7 The Rest- frame Stromgren Colors vs. Red 
Colors 

We have demonstrated in the preceding sections a good cor- 
relation between rest-frame Stromgren colors and param- 
eters determined from the blue part of the spectrum (Hg 
and D n (4000)). In this section we analyze the correlation 
between rest-frame Stromgren colors and parameters deter- 
mined from the red part of spectrum. As shown in Sec- 
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Figure 24. The top panel shows the rest-frame r — i color as a 
function of redshift. The upper middle panel shows the correlation 
between the rest- frame r—i color and the PI rest-frame Stromgren 
color. The dashed line is a best-fit third order polynomial (see 
text). The distribution of residuals between the fit and measured 
values is shown in the lower middle panel by the thick dashed line. 
The thin solid line is a Gaussian with the width of 0.048 mag. 
The bottom panel shows the residuals as a function of the P2 
rest-frame Stromgren color. 
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tion 13.31 the age-metallicity degeneracy cannot be broken 
using rest-frame Stromgren colors employed here. It is pos- 
sible that additional information contained at the redder 
wavelengths could provide a method to separate the age and 
metallicity effects. If the properties of the red part of the 
spectrum can be fully predicted using rest-frame Stromgren 
colors, then such a possibility can be ruled out. Furthermore, 
if the blue and red parts of the spectrum are dominated by 
uncorrelated stellar populations, it should not be possible to 
accurately predict r — i color using blue Stromgren colors. 

To characterize the red part of the spectrum, we syn- 
thesize the rest-frame r — i color, using the same method as 
for synthesizing rest-frame Stromgren colors. The top panel 
in Fig. 1241 shows the rest-frame r — i color as a function 
of redshift. Reassuringly, there is no strong dependence on 
redshift - the upper limit on the change of the median color 
with redshift is ~0.02 mag (in the 0.02-0.12 redshift range). 
The correlation between the rest-frame r — i color and the 
PI color is shown in the upper middle panel. The dashed 
line is a best- fit 

r — i — 0. 109 PI 3 - 0.221 PI 2 + 0.288 PI + 0.238, (12) 

which predicts the rest-frame r — i color with an rms of 
only ~0.05 mag. The distribution of residuals is shown in 
the lower middle panel, and is well described by a Gaussian 
distribution. 

Since P2 is a measure of dust reddening in a galaxy, 
some correlation between the r — i residuals and P2 could 
be expected, although the errors in the determination of PI, 
P2 and r — i, which are comparable to the magnitude of r — i 
residuals, may easily mask it. The dependence of residuals 
on the P2 color, shown in the bottom panel, gives a hint of 
such a correlation. For example, the median r — i residual is 
0.03 mag for P2=0.05, while it is -0.02 mag for P2=-0.05. 

This one-dimensional parameterization of galaxy spec- 
tral energy distributions (SEDs) extends to even longer 
wavelengths. Obric et al. (2004, 2006) find that SDSS imag- 
ing magnitudes in the u and r bands, together with redshift 
information, can be used to predict 2MASS K band (2.2 fim) 
magnitudes with an rms scatter of only ~0.2 mag, with a 
significant fraction of this scatter contributed by measure- 
ment errors (they estimate that astrophysical scatter is much 
smaller than 0.1 mag). In summary, the UV to IR SEDs of 
galaxies can be described as a single-parameter family with 
an accuracy of 0. 1 mag, or better; in addition, the more de- 
tailed spectral diagnostic parameters, such as line strengths, 
appear well correlated with the overall SED. 



5 THE CORRELATION BETWEEN THE 

REST-FRAME STROMGREN COLORS AND 
SPECTRAL EIGENCOEFFICIENTS 

In the previous Section we discussed correlations between 
the rest-frame Stromgren colors and various other spectral 
parameters. Here we study the correspondence between a 
classification scheme based on spectral eigencoefficients pro- 
posed by Yip et al. (2004) and the principal colors PI and 
P2. 

Yip et al. used the Karhunen-Loeve transform to clas- 
sify 170,000 SDSS galaxy spectra. These spectra are classi- 
fied in a plane spanned by the mixing angles between the 
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Figure 25. The bottom panel shows a spectral classification 
plane for galaxies proposed by Yip et al. (2004). Each galaxy is 
represented by a dot color-coded according to its position in the 
P1-P2 plane, shown in the top panel. The sequence from early- 
type to extreme emission line galaxies adopted from Yip et al. 
(2004) is shown on the right side in the bottom panel. Note the 
good correspondence between PI and P2 colors and the Yip et 
al. classification. 



first three eigencoefficients, &kl and Okl- This classifi- 
cation discriminates between early-, late-type and extreme 
emission line galaxies (see their Section 5 for details). Given 
correlations between various spectral parameters discussed 
here, it is expected that the distribution of galaxies in the 
Qkl-Qkl plane is, at least to some extent, correlated with 
the principal colors PI and P2. 

Fig. l25l shows this distribution, color coded according to 
the position of each galaxy in the P1-P2 plane. The corre- 
lation between the PI color and the proposed classification 
is evident. As discussed in detail by Yip et ah, the mean 
spectra for galaxies selected from small ($a'l, Qkl) regions 
correspond to different spectral types. This suggests that the 
PI color is a robust measure of galaxy type, and confirms 
that most of the variance in galaxy SEDs is captured by a 
small number of eigencomponents, as found by Yip et al. 

The correlation with P2 is not as strong as for PI, but 
appears discernible for early- to late-type galaxies. There is a 
slight gradient from purple (P2 < 0) to red-orange (P2 > 0) 
in the 80 < Qkl < 100 range. 
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6 THE REST-FRAME UV COLORS OF z > 0.18 
GALAXIES 

Our analysis of the rest-frame Stromgren colors in the pre- 
ceeding sections did not utilize the uz magnitudes. Since the 
SDSS spectra cover only wavelengths longer than 3800 A, 
the uz magnitudes can be synthesized only for galaxies at 
redshifts greater than z ~ 0.18 (hereafter "high-z" sample). 
The "high-z" sample includes only ~ 9% of the 99,088 galax- 
ies from the "main" galaxy sample, and is strongly biased 
towards red galaxies because of the severe luminosity cut 
(see Fig. \2\\ . Nevertheless, it includes a sufficient number 
of blue galaxies to demonstrate a correlation between the 
uz — vz color and the PI color, shown in Fig. 1261 
The best-fit relation shown in Fig. 1261 

uz — vz = 0.60 PI + 0.77, (13) 

can be used to "predict" the uz — vz color with a resid- 
ual scatter of only 0.08 mag! It is noteworthy that the 
uz — vz residuals are not correlated with the P2 color. 
This small scatter means that either the dependence of the 
uz — vz color on metallicity is not as strong as assumed a 
priori, or that age and metallicity for galaxies are combined 
in such a degenerate way that the distribution of galaxies 
in the uz — vz vs. PI plane is essentially one-dimensional 
(within this 0.08 mag scatter), similarly to their distribu- 
tion in the PI vs. P2 plane. 

In order to distinguish between these two possibilities, 
we resort again to the Bruzual & Chariot population synthe- 
sis models (see Section ^, 31 . Fig. l27l shows model predictions 
for the dependence of the uz — vz color on age and metallic- 
ity. At least for the smallest (<30 Myr) and largest (>2 Gyr) 
ages, the uz — vz color appears sensitive to metallicity. How- 
ever, this "high-z" sample doesn't include extremely blue 
galaxies with implied ages below 30 Myr, or Pl<-0.5 (see 
Fig. |HJ , so model predictions are not at odds with the corre- 
lation displayed in Fig. 1261 The situation is less clear for old 
red (>2 Gyr, or Pl>0.2) galaxies. Model tracks for 10 Gyr 
old galaxies (Pl>0.4) predict a >0.2-0.3 mag wide range 
of uz — vz colors, while the observed distribution has the 
root-mean-square scatter of 0.08 mag. 

The resolution of this discrepancy is the degenerate de- 
pendence of the PI and uz — vz colors on age and metal- 
licity. For galaxies older than ~2 Gyr, both colors become 
redder as age and metallicity increase, by about the same 
amount. As shown in Fig. 1281 for Pl>0.4 the model tracks 
converge and become a single track irrespective of metallic- 
ity. Taken at face value, this convergence suggests that the 
uz — vz color scatter of elliptical galaxies (at a given PI) 
is due to age spread, and not due to a finite width of their 
metallicity distribution. 

This interpretation has to be taken with caution be- 
cause models never produce sufficiently red uz — vz colors 
to explain the observations of the reddest galaxies. While 
one is tempted to invoke the dust (which is not incorporated 
in these models) as an explanation of this discrepancy, these 
galaxies are the reddest and most luminous galaxies which 
are not supposed to have significant amounts of dust. Fur- 
thermore, the models show significantly different behavior 
for blue galaxies than displayed by the data distribution. 
It is quite plausible that model spectra in the rest-frame 
UV wavelength range require significant improvements. For 
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Figure 26. The uz — vz color vs. PI color-color diagram for 
galaxies (dots) with synthesized uz magnitudes (0.18 <z<0.23). 
Large circles show medians for PI bins. The middle (straight) 
dashed line is the best-fit to these medians, and the other two 
dashed lines mark the 2<r envelope around the medians. The best- 
fit relation can be used to "predict" the uz — vz color from PI 
with a residual scatter of only 0.08 mag. 
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Figure 27. The Bruzual & Chariot model predictions for the de- 
pendence of the uz — vz color on age and metallicity ([Z]=0.004, 
squares; 0.02, circles; 0.05, dots) of single stellar populations. Note 
that for small ages (<30 Myr) metallicity can have a strong im- 
pact on the uz — vz color, while for intermediate ages (30 Myr 
<age < 2 Gyr) it plays only a minor role. For old stellar popu- 
lations (>2 Gyr), the uz — vz color becomes again increasingly 
sensitive to metallicity. However, not even for the oldest galaxies 
the model uz — vz color is sufficiently red to explain the observa- 
tions. 



example, Fig. 1281 indicates that a model spectrum fit to 
A > 4000 A spectral range for a blue galaxy (say, Pl=— 0.5), 
will have the flux in the 3200-4000 A wavelength range un- 
derestimated by ~50% (0.5 mag). Note that this discrep- 
ancy cannot be easily attributed to dust effects (which are 
not accounted for in model spectra) because model spectra 
are redder, and not bluer, than the observed spectra. 
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Figure 28. The uz — vz color vs. PI color-color diagram. The 
observed galaxy distribution is shown by linearly spaced con- 
tours (the same distribution is shown by dots in Fig. 1261 . The 
Bruzual & Chariot model predictions are shown by lines and sym- 
bols, analogously to Fig. 1271 The position along each model track 
is controlled by the age of a single stellar population, while dif- 
ferent tracks correspond to different metallicities. Note that the 
effects of age and metallicity are largely degenerate such that dif- 
ferent model tracks do not deviate from each other by more than 
~0.1 mag. Note also the discrepancy between the model tracks 
and the observed distribution. 



7 DISCUSSION AND CONCLUSIONS 

We analyzed the synthesized rest-frame Stromgren colors 
for ~100,000 SDSS galaxies. This narrow-band photomet- 
ric system is well suited for studying galaxy colors in the 
4000-5800 A wavelength range. Galaxies form a narrow lo- 
cus in the rest-frame Stromgren bz — yz vs. vz — yz color- 
color diagram, with a width of only 0.03 mag. This small 
width, in addition to demonstrating that the errors in syn- 
thesized colors are small, shows that the slope of galaxy 
spectra in the 3200-5800 A wavelength range is practically 
a one-parameter family. This conclusion is independent of 
the details of the employed photometric system. 

Using the Bruzual & Chariot population synthesis mod- 
els, we find that the dependence on age and metallicity 
is fully absorbed in the first principal color axis, PI. This 
color is essentially a measure of the Hubble type (in sense 
of early vs. late type galaxies, rather than detailed mor- 
phological classification). It is not possible to break the age- 
metallicity degeneracy with any of the observables discussed 
here. This may not be possible even when using standard 
methods such as Lick indices, at least in the case of ellip- 
tical galaxies (Eisenstein et al. 2003). As demonstrated by 
Eisenstein et al., the strengths of numerous metal lines, such 
as Mg b, can be accurately predicted from the strength of 
the Hs line (see their Fig. 5); hence, no additional informa- 
tion can be extracted from the measurements of those lines. 
On the other hand, Tremonti et al. (2004) report a very 
tight mass-metallicity relation (0.1 dex scatter in metallic- 
ity at a given mass) for ~50,000 SDSS star-forming galaxies. 
Thus, it may be that age and metallicity are strongly cor- 
related, too, which then would provide a trivial solution to 
the problem of breaking the age-metallicity degeneracy. In- 
deed, such a relation has been claimed for disk stars in our 



Galaxy (Twarog 1980; however, see also Feltzing, Holmberg 
& Hurley 2001). 

The second principal color axis, P2, is independent 
of stellar age and metallicity, at least according to the 
Bruzual & Chariot models. We stress that the definition 
of this color is purely observational, and the models explain 
both its value and narrow distribution width. The P2 dis- 
tribution width (0.03 mag) measured using SDSS data is 
smaller than typical observational errors in rest-frame col- 
ors using other data available to date (>0.05 mag). The P2 
color appears most sensitive to the galaxy dust content, as 
demonstrated by the analysis of model-dependent estimates 
of the dust content by Kauffmann et al. (2003a), and also 
by the redder P2 colors for galaxies detected by the IRAS 
survey. Of course, P2 could be indirectly related to metal- 
licity, if the dust content depends on metallicity (that is, 
the observed P2 color could be correlated with metallicity 
as a result of radiative transfer through interstellar dust, 
rather than due to the dependence of the source function on 
metallicity). 

The principal colors are well correlated with numerous 
other spectral parameters determined from bluer and redder 
spectral ranges, as well as with the position of a galaxy in 
the BPT diagram constructed with emission-line strength 
ratios. These good correlations suggest that not only are the 
overall UV-IR SEDs a nearly one-parameter family, but that 
the same parameter also controls the properties of emission 
lines. Most notably, there is a good correspondence between 
the stellar mass-to-light ratio estimates by Kauffmann et al. 
(2003a) and the principal color PI. 

The low dimensionality of galaxy spectra implied by 
the numerous correlations among the continuum and spec- 
tral parameters shown here is in good agreement with the 
conclusions by Yip et al. (2004), who applied the princi- 
pal component analysis to SDSS galaxy spectra. The same 
conclusion was reached in a detailed study of elliptical galax- 
ies by Eisenstein et al. (2003). These conclusions are valid 
for the overwhelming majority of all galaxies; however, we 
emphasize that the analysis presented here is based on low- 
order statistical moments (such as root-mean-square scatter) 
and cannot rule out a possibility that as many as 5-10% of 
galaxies behave differently. 

Our results support the suggestion by Rakos and col- 
laborators (Fiala, Rakos & Stockton 1986) that the rest- 
frame Stromgren photometry provides an efficient tool for 
studying faint cluster galaxies and low surface brightness 
objects, without the need for time-consuming spectroscopy. 
For example, the stellar mass-to-light ratio can be deter- 
mined from PI with an rms scatter of only 0.10. Since the 
brightest cluster galaxies can be selected from SDSS data 
to redshifts of ~0.8 (e.g. Ivezic et al. 2002), this method 
can be used to quantify the rise in star formation activity 
with cosmic epoch, and the Butcher-Oemler effect, with an 
unprecedented accuracy. We are currently undertaking such 
an observational program. 
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